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Abstract 

Assuming a resonant origin of the quasiperiodic oscillations observed in the X- 
ray neutron star binary systems, we apply a genetic algorithm method for selection 
of neutron star models. It was suggested that pairs of kilo-hertz peaks in the X-ray 
Fourier power density spectra of some neutron stars reflect a non-linear resonance 
between two modes of accretion disk oscillations. We investigate this concept for 
a specific neutron star source. Each neutron star model is characterized by the 
equation of state (EOS), rotation frequency and mass M. These determine the 
spacetime structure governing geodesic motion and position dependent radial and 
vertical epicyclic oscillations related to the stable circular geodesies. Particular 
kinds of resonances (KR) between the epicyclic frequencies, or the frequencies 
derived from them, can take place at special positions assigned ambiguously to 
the spacetime structure. The pairs of resonant eigenfrequencies relevant to those 
positions are therefore fully given by KR,A/, fl, EOS and can be compared to 
the observationally determined pairs of eigenfrequencies in order to eliminate the 
unsatisfactory sets (KR,Af, Q,, EOS). For the elimination we use the advanced 
genetic algorithm. Genetic algorithm comes out from the method of natural se- 
lection when subjects with the best adaptation to assigned conditions have most 
chances to survive. The chosen genetic algorithm with sexual reproduction con- 
tains one chromosome with restricted lifetime, uniform crossing and genes of type 
3/3/5. For encryption of physical description (KR,Af, Q, EOS) into chromosome 
we used Gray code. As a fitness function we use correspondence between the 
observed and calculated pairs of eigenfrequencies. 



1 Introduction 

Recently developed observational techniques provide good quality data from observa- 
tions of quasiperiodic oscillations (QPOs) in black hole and neutron star sources [ 1 1. 
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It is shown that in the case of some neutron star atoll sources (e.g. 4U1636-53 J2)) 
the data can be well fitted by the so called multiresonant total precession model (TJ. 
The fits give high precision values of the neutron star spacetime parameters; in addi- 
tion, in some cases the rotation frequency of the neutron star is measured determined 
almost exactly from the QPO independent measurements. This enables us to put some 
limits on equations of state (EOS) describing the neutron star interior. Here we focus 
our attention of the EOS given by the Skyrmion interaction that are very well tuned to 
the data given by the nuclear physics (3). Using the genetic algorithm method, which 
appears to be very fast and efficient, we select the acceptable EOS from 27 types of the 
Skyrmion EOS selected by other methods 1 3 1, concentrating on the source 4U 1636-53. 
For the other five atoll sources, the method gives similar results. 



The results of recent studies of neutron star QPOs indicate that for a given source the 
upper and lower QPO frequency can be traced through the whole range of observed 
frequencies but the probability to detect both QPOs simultaneously increases when 
the frequency ratio is close to ratio of small natural numbers (namely 3/2, 4/3, 5/4 
in the case of atoll sources studies recently, see [ 1 1). Therefore, the multi-resonant 
orbital model based on the oscillation with Keplerian (z/k) and epicyclic vertical (vg) 
or radial {v r ) frequencies was used to explain the observed data [ 1 1. They are calculated 
assuming the spacetime given by the Hartle-Thorne metric 0[5] . The fitting procedure 
have shown that the best results are obtained using the total precession model, where in 
all the sources the upper frequency v u = and the lower frequency v\ = = ve — v r 
[ 1 ]. We concentrate here on the case of 4U 1636-53, when the mass and dimensionless 
spin of the neutron star are fitted to values (x 2 + 1 precision) in the range 



This interval of allowed values of M and j will be used to test the EOS for the neutron 
star in the 4U 1636-53 source. 

3 The neutron star structure 

In neutron star, the strong gravity (i.e. Einsteins gravitational equations) must be rele- 
vant being crucial for the structure equations. The spacetime geometry is assumed to 
be stationary and axisymmetric, the perturbative approach of Hartle and Thorne [4 1 is 
used. One start with the spherically symmetric space time, given by Schwarzschild- 
like metric 



2 Fitting the QPO data 



M = 1.77 ±0.07, j 
M = 1.84 ±0.07, j 



0.051 ± 0.044, 
0.101 ±0.044. 



(1) 
(2) 



ds 2 



+ r 



-e v dt 2 + [1 - 2m(r)/Vp 1 dr 2 
2 (d<9 2 + sin 2 0{d<j, - Ludt 2 )] . 
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To calculate structure of nonrotating, unperturbed star one has to integrate TOV equa 
tion 

dP _ Gm(r)p (1 + P/pc 2 ) [1 + 47rr 3 P/m(r)c 2 ] 



dr r 2 1 — 2Gm(r)/rc 2 



(4) 



where 



m(r) = / Aur pdr (5) 
o 

is the mass inside radius r. The integration is done outward from the center (for 
given values of central energy density p c ) up to surface (where pressure vanishes). 
One obtains the global properties of non-rotating neutron star as its mass M, radius R 
and the internal characteristic profiles as metric coefficients, pressure, energy density 
and number density of baryons expressed as functions of the radial coordinate distance 
from center. 

The rotational effect is, in the linear approximation, given by the Hartle-Thorne 
metric |4j- 



ds 2 = -e" [1 + 2[h + h 2 P 2 )} At 2 

[l + 2(m + m 2 P 2 )/(r -2M)} 
+ 1 - 2M/r 

+ r 2 [1 + 2{v 2 - h 2 )P 2 ] 

32 i „;„2 



dr 2 



x [d6 2 + sin 2 0(d<f> - ujdt 2 )] , (6) 

where P 2 = P 2 (cos9) = (3cos 2 # — l)/2 is the Legendre polynomial of 2nd order u> 
is the angular velocity of local inertial frame, which is related to star's angular velocity 
£1 and ho,h 2 ,mo,m 2 are functions of r and are all proportional to £! 2 .The angular 
velocity w can be found by solving equation 

1 d ( 4 . dQ \ t 4 dj r 



r 4 dr y ^ dr ) ~^ r dr ~ ^ ^ 

where 

j(r) = e -^ r ^ 2 [l - 2m{r)/r] 1 / 2 . (8) 

One integrate equation (O outward from center for arbitrarily chosen u) c with boundary 
condition dco/dr = 0. At the surface one can calculate the angular momentum J and 
frequency of rotation ft corresponding to ui c from relations 



1 „ 4 fdu 



6 V dr 



J = ^ 3- (9) 



r=R 



2 I 

n ncw = q(r) + —. (io) 

We take frequency of rotation ft as an input parameter, thus after integration of eq. (O 
we rescale the ui to obtain the right value of ft 

w n ew (r) = w old (r) . (11) 

"old 
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After this one should calculate mass and pressure perturbation factors mo and po from 
equations 



dran . 2 di? , _ 1-24' do) 



^ — (E + P)p + -j^[ — 



dr dP K ' r 12 J V dr 

dr 



2 



n r "J" W ! ( 12 ) 



dp m (H-87rr 2 P) 4ir(E + P)r 2 

dr (r — 2m) 2 r — 2m ^° 

1 f 4 j 2 ( duj\ 2 Id/ r 3 j 2 u; 2 



12 r — 2m \_ dr J 3 dr \r — 2m/ 

The mass of rotational object is then given by 

M{R) = M (R) + m (R) + J 2 /R 3 . (14) 

In the next approximation, the quadrupole moment q of the star is introduced. The 
fitting procedure show that q ~ j 2 [1| where j — J/M 2 . Therefore, the spacetime 
could be considered quasi-Kerr and the Kerr metric and related formula of the orbital 
motion could be used. 



4 Skyrmion interactions and realated EOS 

The effective Skyrmion interaction implies a variety of parametrization in the frame- 
work of mean-field theory. All give similar agreement with experimentally established 
nuclear ground states at the saturation density no, but they imply varying behaviour of 
both symmetric and asymmetric nuclear matter when density grows (up to 3no). 

The general form of effective Skyrme interaction implies total binding energy of a 
nuclei as the integral of an energy density functional H, determined as a function of 
nine empirical parameters to, ti, t^, t%, Xo, x\ , X2, X3 and a in the form 

TC = JC + Ho + Ti3 + H-cS, (15) 

where the kinetic term JC — (h/2m) r, is given by the kinetic densities r = t„ + r p , 
with t being the total density, r n (r p )is the neutron (proton) density. The other terms 
are 

given by the relations 

Ho = jt [(2 + x )n 2 - (2x + l)(n 2 p + n 2 n )] (16) 

H 3 = ^t 3 n a [(2 + x 3 )n 2 - (2x 3 + l)(n 2 p +n 2 n )] (17) 
n eff = 18[ti(2 + a; 1 )+t2(2 + a;2)]r„ 

+ - [t 2 (2x 2 + 1) - h (2xi + 1)] (t p Tip -\- T n n n 

) (18) 

8 
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The pressure is then given by 



P(n b ,I) 




2 9e 



(19) 



where e is the binding energy per particle and / = (N — Z)/A denotes assymetry of 
nuclear matter. 

In [3], 87 different Skyrme parametrization were limited to 27, using limits implied 
by the spherically symmetric models of neutron star models and by experimentally 
tested properties of nuclear matter. Here, we shall test acceptability of nine of these 27 
Skyrme parametrization to the limits put by QPO measurements and using axisymmet- 
ric models in first approximations with respect to the star rotation. 

5 Determination of neutron star structure using genetic 
algorithm 

Genetic algorithm (GA) appears from method of natural selection, when subjects with 
best adaptation to assigned conditions have highest chance to survive [6|. GA takes 
into account the following natural mechanism - mutation and lifetime limit restricting 
risk of degradation, which is kept in local extreme from optimization viewpoint. GA 
has iteration character. GA doesn't work with separate result in particular iterations, 
but with population. In each iteration GA works with several (generally a lot of re- 
sults, standard value is hundreds) results, which are included in the population trying 
to ensure appearance of still better results via genetic operations with these results. 
Generally, the GA scheme is given in the form 



where P is population containing N elements, 8 is parent selection operator which 
selects u elements from P. Evaluation for each chromosome performed by the fitness 
function /. 



Genetic operators included in il are namely crossover operator flc, mutation operator 
r^M an d others problem-oriented or implementation-oriented specific operators, which 
all together generate v offspring from u parents. ^ is deletion operator, which removes 
v selected elements from actual population P(t). v elements is add to new population 
P(t + 1) after it, r is stop-criterion. Parent selection operator 6 and genetic operators 
O have stochastic character, deletion operator is generally deterministic. 

We have selected GA with sexual reproduction containing one chromosome with 
restricted lifetime parameter, 5 iterations. The crossover operator is selected with uni- 
form crossing, using genes of type 3/3/5. We used the Gray code for encrypt parameters 
to chromosomes, which is useful by reason of bypassing so-called Hamming barrier. 

Chromosomes are compound of genes; each gene presents 1 bit value. However, 
genes contain more than 1 bit (using the redundancy encrypt). Bit values inside a 
gene are mapped onto outside value of gene (0 or 1) via specific map function, in 



GA = (N, P, f, 0, fi, r) 



(20) 



f:Si*-R, i = l,...N. 



(21) 
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Table 1: Fitness function and corresponding parameters for all tested EOS, £1 € 

(50; 8000) 



EOS type 


Pc 


n 


Fitness 


0-SkT5 


2.1244 


964.17 


56.556 


1-SkO' 


1.5575 


979.70 


56.686 


2-SkO 


1.3987 


711.85 


56.601 


3-SLy4 


1.4896 


1373.71 


57.292 


4-GS 


1.1225 


743.89 


57.223 


5-SkI2 


1.1185 


1241.72 


57.688 


6-SkI5 


0.9359 


791.43 


57.059 


7-SGI 


1.0170 


986.49 


57.028 


8-SV 


0.7889 


590.55 


56.393 



which border between and 1 is not crisp, but there exists so-called "shade zone", 
where value carried by the gene is determined randomly. |J8] [9) .We pressed number of 
members in each generation to 400 (see Fig.l) and interate 40 generations. We define 
two task of selection neutron star structure. The first one determines p c , il with preset 
EOS where fitness function f is \ 2 DO- It is important to say here, that at present the 
fitness function is tabled. The second one determines the EOS, p c , ft with using fitness 
function x 2 - Both tasks put limits on the EOS under consideration. Overcoming those 
limits implies removing the corresponding chromosomes by setting value of fitness 
function to a maximum value (10 200 ). For a given EOS allowed neutron star structure is 
provided the GA described above and fitness function \ 2 - The chromosome structure is 
given by two values, central density and rotation frequency. Interval of central density 
we choose from 0.4 • 10 15 to 3.3 • 10 15 g.cm~ 3 . 

In the first task, partition of central density is set to 4096 (12 bits). Interval of 
rotation frequency is set from 50 to 8000 and partition to 8192 (13 bits). Following 
table shows determination of central density (in units of 10 15 g.crn" 3 ) and rotation 
frequency for given number of EOS. Because the resulting frequencies do not fit the 
observed rotation frequency interval [7|, and also the fitness function has too many of 
local minima with same value for different combination of pc , f2 ,we make determina- 
tion of values mentioned before, but we set range of rotation frequencies from 1827.21 
to 1828.03 with partition 64 (6 bit). 

Clearly, for selecting the most convenient EOS we need referenced % 2 value as 
global minimum x 2 over all EOS. Thus we make second task, determination of the 
best resulting neutron star structure over all used EOS. The corresponding EOS value 
is given by zero based EOS number (4 bits). 

6 Conclusions 

Analysis of QPO's in neutron star atoll sources in the framework of Hartle-Thorne ge- 
ometry gives a very detailed data (neutron star parameters as mass, spin and quadrupole 
momentum) that could be quite well used for constraining the wide scale of allowed 
EOS constraining the structure of neutron stars by solving the complex set of structure 
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Table 2: Fitness function and corresponding parameters for all tested EOS, £1 G 
(1827.21; 1828.03) 



EOS type 


Pc 


Q. 


Fitness 


0-SkT5 


2.9978 


1827.928 


62.979 


1-SkO' 


1.7896 


1827.453 


57.542 


2-SkO 


1.6762 


1827.210 


57.899 


3-SLy4 


1.5747 


1827.300 


57.611 


4-GS 


1.3084 


1827.236 


57.830 


5-SkI2 


1.2174 


1827.505 


58.061 


6-SkI5 


1.0803 


1827.377 


58.277 


7-SGI 


1.1146 


1827.492 


58.658 


8-SV 


0.8858 


1827.428 


58.45 



Table 3: Global minima for Q e (1827.21; 1828.03) 



^min 




Fitness 


EOS 


Pc 


n 


1827.21 


1828.03 


57.542 


1 


1.79021 


1827.66 



differential equations. The testing by standard approaches is a long time consuming 
procedure. Here, we show in the case of atoll source 4U 1636-53 that the genetic algo- 
rithm method could make the proper selection in a wide sample of EOS of Skyrmion 
type in a very efficient and short way (with same precision and the time consumed for 
the GA about 30 minutes, being by orders shorter that the time consumed by the stan- 
dard methods). We can conclude that using the GA with chromosome (EOS,p c , £1) and 
putting physically motivated restrictions on the angular velocity (£1) and central den- 
sity (p c ) governing the neutron star model being an output of the structure equations 
and the fitness function x 2 m the algorithm procedure, we are able to find the most 
probably structure of neutron star that fits the observed QPO data, with respect to ob- 
served rotational frequency for the source 4U 1636-53. It is given by the chromosome 
(1, 1.79- 10 15 , 1837.66). 

We expect that the method allows much strongest test by the genetic algorithm 
including the quadrupole momentum calculated directly without using the quasi-Kerr 
approximation q ~ j 2 . 
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